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Abstract 

In this paper, we present a new approach to constructing a “self-avoiding” walk t ‘ iro ’^ 
a triangular mesh. Unlike the popular approach of visiting mesh elements using spa fi g 
based on a geometric embedding, our 

that it uses the mesh connectivity only. We present ar .algorithm f« -onnottMMW » 
avoiding walk which can be applied to any unstructured triangular mesh. The complexi y 
is 0(» • log(n)), where n is the number of triangles in the -h. We^w 
that for hierarchical adaptive meshes, the algorithm can be easily parallelized by t g 
advantage of the regularity of the refinement rules. The proposed approach should be very 
useful in the runtime partitioning and load balancing of adaptive unstructured grids. 


1 Introduction 

Advances in adaptive software and methodology notwithstanding, parallel computational strate- 
git „Tbe an Len.ial ingredient in solving complex, real-life problems. However paraM 
Lmputers are easily programmed with regular data structures; so the f 

parallel adaptive algorithms for unstructured grids poses a serious challenge An efficient par 
allelization of these unstructured adaptive methods is rather difficult, P n man y ne o t le 
imbalance created by the dynamically-changing nonumform and .rrege dar g ^ ^oneR>e J 
it is generally believed that unstructured adaptive-grid techniques will constitute a sigmhcant 
fractfon of future high-performance supercomputing. Mesh adaption and dynamic load balanc- 
ffig ^r S t be accomplished rapidly and efficiently, so - not to cause a significant overhead to 

the numerical simulation. i 

Serialization techniques play an important role in all parts of a Finite Element Method 
(FF,M) over adaptive unstructured grids. A numbering of the unknowns allows or ■a ^matrix- 
vector notation of the underlying algebraic equations. Special numbering 

McKee frontal methods) have been developed to optimize memory usage and locality of the 
algorithms” On the other hand, in many cases, runtime support for decomposing dynamic 
adaptive grid hierarchies is based on a linear representation of the grid hierarchy [8] in e or 
of a space-filling curve. Salmon et al. [10, 11] have demonstrated the successful appl.eat.on o 
techniques based on space-filling curves to N-body simulations. Other researchers [4, 6, , , ] 

have also shown how space-filling curves can be used for graph partitioning and sum ar grap - 

,e '“he P g r“rndca of a space-filling curve is a sort of a serialization (or linearization) of visiting 
points in a higher-dimensional space. A standard method for the construction ,* to embed 1 the 
object of study into a regular environment (where the standard space-tilling curves live). . 
approach introduces a kind of an -artificial” structure in the sense that the entire construction 
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depends on the embedding. 1 Furthermore, this type of a linear representation forgets about 
the combinatorial structure of the mesh which drives the formulation of operators between the 
finite element spaces. The questions arises whether one could reduce the different requirements 
for a serialization in an adaptive FEM over unstructured grids to a common denominator. 

In this paper, we present a new approach to constructing a “self-avoiding walk' through 
a triangular mesh. Unlike the popular approach of visiting mesh elements using space-filling 
curves which is based on a geometric embedding, our approach is combinatorial in the sense 
that it uses the mesh connectivity only. We present an algorithm for constructing a self-avoiding 
walk which can be applied to any unstructured triangular mesh. The term unstructured refers 
to information that can hardly be compressed, and which does not contain any symmetry or 
redundancy. At this level, any labeling or numbering of the mesh components is as good (or as 
bad!) as any other scheme. However, the situation changes if the mesh adaption is done in a 
hierarchical manner where we, based on a set of simple rules for coarsening and refinement, pay 
particular attention to the triangulation of the initial mesh. In contrast with a general adaption 
scheme that produces just another unstructured mesh, this strategy exploits the structure of 
the adaption hierarchy to simplify the labeling process. We would like to modify (or adapt) the 
labeling only in regions where the mesh has been altered but not rebuild the entire indexing. 
This idea has been extensively developed in [3] . 

Self-avoiding walks can be used to improve the efficiency of the respective algorithms, in 
particular locality and load balancing. However, we need to search for some special c asses 
of self-avoiding walks that will facilitate runtime load balancing with good locality in terms 
of memory access or communication. Since we are not interested in just any kind of walks, 
we have to first tackle the question of the existence of walks with the desired properties. In 
this paper, we give an existence proof which is constructive in nature so that it also provides 
an algorithm for generating such walks. For the case of hierarchical adaptive refinement, the 
proposed algorithm is easily amenable to parallelization. In the process of deriving the existence 
of a particular class of walks, a natural mathematical framework emerges which, we believe, 
can be easily adapted for othei* special classes of walks. 

In Section 2, we define a certain class of walks whose existence for arbitrary unstructured 
meshes we can prove (see Appendix A). In Section 3, we present an algorithm based on the 
existence proof and show some sample results. The algorithm works for arbitrary unstructured 
meshes; however, a parallelization is non-trivial and the walk has to be completely rebuilt after 
mesh adaption. Furthermore, since we cannot make any regularity assumptions, the chances are 
quite low that we could prove the existence of constrained (in the sense of boundary conditions) 
walks. In Section 4, we briefly review hierarchical mesh adaption and a suitable indexing 
scheme. We find the missing regularity that enables us to prove the existence of constrained 
walks (see Appendix B) and to formulate a truly parallel algorithm for the construction of walks 
which is well-behaved with respect to mesh adaption. Section 5 concludes the paper with some 

perspectives. 


> It is similar to describing metric properties of a curved surface in terms of a parametrization in the envel^Pmg 
(uncurved) three space instead of doing inner geometry, i.e., geometry based only on measurements on thesur - 
itself. The only inner properties of a mesh are its decomposition into cells of different dnnens.ons and a face 

relation describing how the cells are glued together. ... , , , r 

2 This term should be familiar to people working in Monte Carlo methods; however, it is beyond the scop . f 

this paper to address the similarities and differences in detail. 
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2 Definitions 


Consider an arbitrary two-dimensional triangular mesh 3 9JT and denote the underlying set of 
triangles by H. Let #H denote the cardinality of H. 

Definition 2.1. A mapping 

W :{1,. (!) 

is called a self-avoiding walk 4=^fief w **' a bijection and 

Vi € {1, • • • , “ 1} w(i) an d + 1) s ^ are an e ^ e or a vertex ■ 

Remark. The previous definition states that a walk visits each triangle exactly once and 
that jumps (i.e., two consecutive triangles in the walk that have an empty intersection) are 
forbidden. The reader may wonder why our requirements for a walk (cf. Eq. (2)) are rather 
weak As simple examples show, the condition that triangles following one another must s are 
an edge, is too strong. Generally, self-avoiding walks do not exist under this assumption (cf. 

Fig. 1 for a trivial example). 



Figure 1: A simple counterexample (look at the dual graph!). 

In the following, we consider a special class of self-avoiding walks. 

Definition 2.2. A self-avoiding walk is called proper •< — >def 

Vi € {2,... ,#R- 1} w(i-l)n4)^w(i)nw(i + l). (3) 

Remark. The previous definition states that for three triangles following one another in a 
proper self-avoiding walk, “jumping” twice over the same vertex is forbidden (cf. Fig. 2). 



Figure 2: A forbidden jump. 

3 By a mesh we understand a simplicial complex coining from the simplicial decomposition of a connected 2D 
manifold (with or without a boundary). 
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3 Existence 

Proposition 3.1. There exists a proper self-avoiding walk for an arbitrary triangular mesh 

©t. □ 

For a proof of Proposition 3.1, we refer the reader to Appendix A. The proof is inductive 
over the number of triangles in the mesh and extends an existing proper self-avoiding walk over 
to a larger mesh. The proof provides a set of elementary rules which can be used to formulate 
an algorithm for constructing proper self-avoiding walks over arbitrary unstructured meshes 
The value of such rules becomes apparent if practical questions such as the proper extension of 
existing “incomplete” walks are addressed (and the difficulty of the proof is in this part). 


3.1 The Basic Algorithm 

The proof of Proposition 3.1 provides an algorithm for the construction of proper self-avoiding 
walks: starting from an arbitrary triangle of 9Jt, choose new triangles sharing an edge with the 
current “subcomplex” and extend the existing proper self-avoiding walk over the new triangle. 
This algorithm can be represented recursively. Figure 3 shows the pseudo-code for the basic 
function called EXTEND-WALK. The algorithm begins with the selection of an arbitrary 
triangle from the mesh. After removing this triangle from the mesh and inserting it (as the first 
triangle) into the walk, the function EXTEND-WALK is called with P being the position of 

the first triangle in the walk.' 1 


EXTEND_WALK(M, W , P) 

mesh M, walk W, position of a triangle in walk P; 

while (triangle T remains in M sharing an edge with the triangle at P) { 

insert T into W according to the rules discussed in the proof of Proposition 3.1; 
remove T from M; 

set new position Q to position of 2 in IT, 

EXTEND -WALK (A/, W , Q); 

} 

} 


Figure 3: The basic function. 


The complexity of this algorithm is 0(n • log(n)), where n is the number of triangles in the 
mesh This complexity is obtained as follows. Since we organize our triangles in a red-black 
tree, a search in this structure is of complexity 0(log(n)). The edge-triangle incidence relation 
can also be organized in a red-black tree. Since we have to perform this search for each triangle 
in our mesh, the overall complexity of our algorithm is 0(n • log(n))- 

Given the existence of proper self-avoiding walks, other approaches might be favorable For 
example, self-avoiding walks have a long tradition in the application of Monte Carlo methods 
to study long-chain polymer molecules. However, this is not the focus of this paper. 


+ u a position HI a walk, we mean a reference to a triangle in the walk rather than an index. We need a data 
JilSi,. b-rt*;, therefore, a.. ore,,, » no. appropriate. I„ C++ „„ple,n.n.ario„. a porerion » 
given by list<Triangle> : : iterator. 
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3.2 Results 

We implemented a sequential algorithm for the generation of a proper self-avoiding walk using 
the set and map containers from the C++ Standard Template Library winch* part of the 
ANSI-C++ standard. Both containers are implemented with red-black trees. r 

presented in Tables 1 and 2 were obtained for some sample meshes on a 110 MHz microbFAKG 
II using GNU’s g++ compiler (version 2. 7.2.3) with the -0 ilag. Note that t le times in e 
tables below include the setup time for the mesh and the edge-triangle incidence structure. 


triangles 

edges 

ratio 

time (s) 

65536 

131073 

0.500 

12.85 

65536 

114690 

0.571 

12.27 

65536 

106500 

0.615 

12.21 

65536 

102408 

0.640 

12.15 

65536 

100368 

0.653 

12.13 

65536 

99360 

0.660 

12.02 

65536 

98880 

0.663 

11.99 

65536 

98688 

0.664 

11.82 


triangles 

edges 

time (s) 

2048 

3136 

0.30 

4096 

6240 

0.57 

8192 

12416 

1.31 

16384 

24768 

2.71 

32768^ 

49408 H 

5.54 

65536 

98688 

11.82 

131072 

197120 

24.07 


Table 1: Runtimes for various mesh sizes. Table 2: Runtimes for various triangle/edge ratios. 


The results in Table 1 indicate the dependence of the runtime on the number of triangles 
(or edges). According to our description of the algorithm, the runtime may also depend on 
the triangle/edge ratio. The results in Table 2 explore this possibility; however, the numbers 
clearly show that the dependence of the execution time on the triangle/edge ratio is negligible 
and confirm our complexity estimate for the whole algorithm. 


4 Optimization and Parallelization 


A standard parallelization of the deterministic algorithm would be based on the observation 
that if the proper extensions of a self-avoiding walk for two triangles do not interfere with 
each other, then they can be done in parallel. Such an approach is feasible in dealing with 
single meshes, when no additional structural information is available, or with families of meshes 
without any hierarchical structure, although this procedure might be computationally expensive 

in connection with adaption. r . r , t 

For a extended discussion of unstructured mesh adaption, we refer the reader to [ , oj. n 
the following, we restrict ourselves to the case of hierarchical refinement and apply the standard 
1:4 (red) and 1:2 (green) subdivision rules for triangles (cf. Fig. 4). 

A couple of additional rules are then applied, primarily to assure that the quality of the 

adapted mesh does not deteriorate drastically with repeated refinement: 

1. All triangles with exactly two bisected edges have their third edge also bisected. Thus, sue i 

triangles are isotropically refined. ... , , 

2. A green triangle cannot be further subdivided. Instead, the previous subdivision is discarded 

and isotropic subdivision is applied to the (red) ancestor triangle. 


5 One might also think about an implementation using hash tables. 




Figure 4: The left picture shows the isotropic subdivision of a triangle. The right one gives an 
example of anisotropic subdivision. 


For hierarchical meshes a completely different approach to constructing proper self-avoiding 
walks is possible: the general algorithm has to be applied only once to the initial mesh and this 
“coarse” self-avoiding walk can be reused after adaption. On the local level, we are going o 
exploit the regularity of the refinement rules. This leads us to the consideration of constrained 
self-avoiding walks. Local regularity allows us to prove the existence of solutions for the un- 
derlying “boundary value problems”. (In that sense, the algorithm der.ved from the proof of 
Proposition 3.1 guarantees the solvability of “initial value problems” only.) In the next section, 
we sketch a technique which allows the efficient handling of constraints. 


4.1 An Indexing Technique for Hierarchical Meshes 

It is the task of an index scheme to properly name or label the various objects (vertices edges, 
triangles) of a mesh. We prefer the term index scheme instead of numbering to stress that the 
use of natural numbers as indices is not sufficient to meet the naming requirements of the FE 
objects on parallel architectures. For a detailed discussion, refer to [2, 3]. 

Our index scheme is a combination of coarse and local schemes. The coarse scheme labels 
the objects of the coarse mesh in such a way that the incidence relations can be easily derived 
from the labels. The vertices are enumerated starting from 1. For the example in 1 ig. o, e 


1 2 



Figure 5: An L-shaped domain and its coarse triangulation, 
set of vertices for the coarse triangulation consists of the following numbers. 


vertices = {1, 2, 3, 4, 5, 6, 7, 8}. 

The edges of the coarse triangulation are indexed by ordered pairs of integers that correspond 
to the endpoints of the edges. The ordering is chosen so that the first index is less than t e 


second one. For the example in Fig. 5, the set of coarse edges consists of the following pairs: 

edges = {(1,2), (1,3), (1,5), (2,5), (3,4), (3,5), (3,6), 

(4,6), (5,6), (5,7), (5,8), (6,7), (7,8)}. 

The same principles are applied to index the coarse triangles. They aie denoted by the trip 
consisting of their vertex numbers in ascending order. Thus, the set of coarse triangles reads: 

triangles = {(1,2,5), (1,3,5), (3,4,6), (3,5,6), (5,6,7), (5,7,8)}. 

Note that this index scheme can be applied to elements with curved boundaries as well. 

The ideal model for local considerations is given by the two-dimensional standard simplex 

a 2 (cf. Fig. 6): 


<t 2 := {(zi, *2,3:3) G M 3 |*i +*2 + 3:3 = I}- 


( 4 ) 


h z 



Figure 6: The two-dimensional standard simplex and some points on it. 

The local scheme exploits the regularity (and the finiteness) of the refinement rules to 
produce names for the objects at subsequent refinement levels [3]. We use (scaled) natural 
coordinates as indices in the local scheme. Again, this is done m a way such that incidence 
relations and the refinement level are encoded in the indices of the objects (cf. ig- )• 
example, the set of vertices at level k in the local model is given by: 

Vk(c j 2 ) := {(a, b,c) £ N 3 | a + b + c = 2*}. ( J ) 

Obviously, there holds v £ V k (a 2 ) ^ 2t» £ and therefore we have the following 

embedding V k {a 2 ) C V k+ i{a 2 ). This shows that we can easily move between refinement levels 

^ Dmmte S by E k (a 2 ) the set of edges at level k. We choose the integertriple corresponding to 
the midpoint of an edge as its index. Hence E k {o 2 ) ^ Vfe+i(cr 2 ) - V k (a ) and 


k 

pj(k) _ (cr 2 ) := | | Ei(& 2 ) =+ ~ Vfc + i(n 2 ) — Vq{<t )■ 

i=o 


( 5 ) 
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( 0 , 0 , 1 ) 


( 0 , 0 , 2 ) 




Figure 7: Examples of the local index scheme for triangular elements. The vertices and triangles 
are denoted by integer triples (triangles by bold face). 


Denote by T k {a 2 ) the set of triangles on level k (red and green) and let 

T (fc) = t W {o 2 ) := \_\Tiia 2 ). ( 7 ) 

i = 0 

We choose the integer triples corresponding to the barycenter of a triangle as its index. 

The coarse and local schemes are combined by taking the union of the Cartesian products 
of the coarse mesh objects with their corresponding local schemes. Ambiguities are resolved by 
using a normal form of the index. The key features of such a scheme are. 

• Each object is assigned a global name that is independent of any architectural considerations 
or implementation choices. 

• Combinatorial information is translated into simple arithmetic. 

• It is well-behaved under (adaptive) refinement. No artificial synchromzation/serialization is 
introduced. 

• It can be extended (with appropriate modifications) to three dimensions [3J. 

4.2 Constrained Self- Avoiding Walks 

We call a subset T C T^(c 7 2 ) a local refinement of <r 2 if it can be obtained by applying the 
refinement rules described in section 4. 6 The refinement level of a coarse triangle is defined 
to be zero. For an arbitrary triangle, the refinement level is defined as the successor of the 
refinement level of the parent triangle (the triangle it was created from by subdivision). The 
refinement level of a triangulation is defined to be maximum of the refinement levels of its 
triangles [3]. We denote the level of a triangulation T by l(T). 

The scheme developed in the previous section allows us, for the case of hierarchical re- 
finement, to decompose our considerations into global and local cases. Here “global means 
“related to the initial (coarse) mesh” , while “local” means “restricted to a triangle of the coarse 
mesh”. Given a walk over and a local refinement of the coarse mesh one might ask whether it is 
possible to extend the walk over the triangles of the local refinement. In addition, we would like 
to decouple the considerations for different coarse triangles. This leads naturally to what we 
call constrained self-avoiding walks: a walk over the initial mesh leaves a trace or footprint on 

6 S i nce wc allow coarsening, the term “local refinement” is somewhat misleading. However any local refinement 
obtained by a combination of coarsening and refinement can be obtained by refinement only, discarding the history 
of the adaption. This assumes that coarsening did not occur beyond the initial mesh. 
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each coarse triangle which is then translated into a “boundary value problem” for the extension 
of the self-avoiding walk. The regularity of the local picture allows us to prove the existence of 
constrained proper self-avoiding walks. 

Definition 4.1. A pair (a, ft) G Vi(cr 2 ) x V\ (a 2 ) is called constraint on a 2 4=4* de f (5. 

Remark: The previous definition states that we forbid constraints requiring that a walk enters 
and leaves a triangle through the same edge or the same vertex. (The latter would imply 
a forbidden jump in the “coarse” walk.) In the definition, we implicitly used the embedding 
V 0 (a 2 ) C Fi(a 2 ) and the isomorphism E 0 {a 2 ) ~ V^-Voia 2 ). Figure 8 shows a few examples 
of constraints. The entry constraints are denoted with bullets and the exit constraints with filled 
triangles. 



Figure 8: A few examples for constraints. 


Definition 4.2. Let c = («,/?) be a constraint on a 2 and T C T^(a 2 ) a local refinement of 
a 2 . A proper self-avoiding walk to : {1,... , #T) -> T is called compatible with c ^ def 

a C oi(l) and j3 C iofifiT). (8) 

Lemma 4.1. Let c = (a, (5) be an arbitrary constraint on a 2 . Let R,G G T\(a 2 ) be a red 
and green triangle, respectively. There exist proper self-avoiding walks on R and G which are 

compatible with c. □ 

Remark: Lemma 4.1 can be easily proved and is valid not only for constraints on a 2 and 
R,G G T l (<t 2 ), but also for constraints on red and green triangles at any level, since the proof 
of Lemma 4.1 depends only on the fact whether a triangle is red or green. 

Obviously, the formulation of constraints is not restricted to coarse triangles, i.e., triangles 
at level zero. A slightly more general definition than definition 4.1 goes as follows: 

Definition 4.3. Let k G N, k > 0. A pair (a, ft) G Vfc(<r 2 ) x Vfc(<7 2 ) is called constrain t on 
t G Tk-\((J 2 ) < c d ef a P- 

Proposition 4.1. Let c = ( a, (3 ) be an arbitrary constraint on a 2 and T C T^ k \o 2 ) be a local 
refinement of a 2 . There exists a proper self-avoiding walk io T , c on T which is compatible with 

c. □ 

Given a constraint on a 2 and a local refinement T C T«(< r 2 ), Proposition 4.1 guarantees 
the existence of a compatible proper self-avoiding walk. The proof in Appendix B also provides 
an algorithm for constructing the solution of this “boundary value problem”. The synthesis 
of Propositions 3.1 and 4.1 provides a simple parallel algorithm which requires changing an 
existing proper self-avoiding walk only where the mesh has been adapted (coarsened/refined). 
If the mesh adaption affected only levels higher than k, the walk has to be modified only over 
triangles whose level is greater than k. Note that the local adaptions of the walk are decoupled 
from one another and can be done independently; thus, this is a truly parallel algorithm. 
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5 Conclusions and Perspectives 

In this paper, we have developed a theoretical basis for the study of self-avoiding walks over 
two-dimensional adaptive unstructured grids. We described an 0(n • log(n)) algorithm for the 
construction of proper self-avoiding walks over arbitrary unstructured meshes and reported 
results of a sequential implementation. We discussed parallelization issues and suggested a 
significant improvement for hierarchical adaptive unstructured grids. For this situation, we 
have also proved the existence of constrained proper self-avoiding walks. 

The algorithms presented allow a straightforward generalization to tetrahedral meshes; how- 
ever, there is some flexibility because of additional freedom in three-space. The application to 
load balancing is likewise rather straightforward. However, modifications to improve locality m 
the sense of a generalized frontal method requires a more detailed study of the relationships 
between (more specialized) constraints on walks and cache memory behavior. 
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A Proof of Proposition 3.1 

Before we start with the proof of Proposition 3.1, we remind the reader of two technical results. 
The respective proofs axe not very difficult, but beyond the scope of this papet. 

Lemma A.l. Let t,t u t 2 ,t 3 be triangles in a mesh 911 and 


Vi € {1,2,3} tj|i. 

(9) 

Then either 


ti n t 2 n t 3 = 0 

(10) 

or DJI is tetrahedral (cf. Fig. 9). □ 



(Note that if two triangles in a mesh share two vertices they must share the corresponding 
edge.) 



Figure 9: A tetrahedral mesh. 


Lemma A.2. Let 911 be a mesh. There exists a triangle r in 9tt such that 911 t, the complex 
obtained by removing r from 911, is still a mesh. □ 

We use the following notation. If two triangles t u t 2 share an edge, we write t x \ t 2 . Let w 
be a self-avoiding walk. If a >(i) | w(i + 1), we will write w(i) h w(i + 1) indicating that w enters 
w (i + 1) from a ;(t) over an edge. If w(i) \ w(i + 1), we write w(i) rv w(t + 1) indicating that w 
jumps into u>(i + 1) from u)(i) over a vertex. 

Proof of Proposition 3.1: We prove Proposition 3.1 by induction over the number of trian- 

glef> Assume that there exist proper self-avoiding walks for meshes with n triangles. Let 911 be 

a mesh with n + 1 triangles. „ , . . . rtW 

Let r be a triangle in 911 such that 911 - r, the mesh consisting of all the triangles in 911 

except r, is a mesh (see Lemma A.2). Mesh 911 - r has n triangles and hence, by our induction 
assumption, there exists a proper self-avoiding walk wot-t- We show that ujm-r can be exten e 
to a proper self-avoiding walk u*m for 911. We call this a** a proper extension of wan-r* 

Let t be an triangle of 9H - r sharing an edge with r (i.e. t \ t holds). For an appropriate 

i £ { 1, . • • , n}, there holds 

t = WOT -A*)- 

The discussion naturally can be split up into four cases, depending on whether wot -t enters 
t. through an edge or a vertex and leaves through an edge or a vertex, respectively. 
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We omit the subscript - r from u> in the following for the sake or clarity. If we do not 
know whether the transition from w(i) to w(i + 1) goes over an edge or a vertex, we write 

^ In tlieVigures below, parts of o*»i_ T are drawn as solid lines whereas the modifications leading 
to are drawn as dashed lines. Triangles are indexed by their position in the walk (e.g. w(») 
is denoted by i ) . 


Case I: u( i — 1) b v{i) b + 1) 

Figure 10 illustrates the modifications necessary if enters and leaves t throng an 

edge. 




(b) 



Figure 10: Existing walk enters the triangle adjacent to r through an edge and leaves it through 
another edge. 


• If u >(* - 1) n r ^ u>(i — 2) n u){i — 1), then u>(i — 1) -> r h uj(i) is a proper extension of u 
(cf. Fig. 10(a)). 

• If w(* - 1) n r = w(i - 2) n w(t - 1), then a >(i -2) rx u>(i - 1). 

- If ui(i + 1) fir / a >{i + 1) fl u >(i + 2), then ai(i) hr -4 u>(i + 1) is a proper extension 
of u> (cf. Fig. 10(b)). 

- If uj{i + 1) D r = u>(i + 1) n u>(i + 2), then w(i + 1) rv u){i + 2). 

* Ifw(*-l)nu;(t + l)^w(t + l)nta»(i + 2),thfino;(t-2) -+ r F w(t) H w(* - 1) 
u){i + 1) is a proper extension of u> (cf. Fig. 10(c)). 

* Let w(i - 1) n u >(i + 1) = w(t + 1) n w(t + 2). According to Lemma A.l, either 
SOT is tetrahedral or w(t - 1) fl w(t + 1) D r = 0- If 3H is tetrahedral, there is 
nothing to prove. For the other case, we have 0 = u>(i — 1) H u>(i -(- 1) fl r — 
w (i + 1) n w(« + 2) Hr = w(t + 1) n r = w(i + l)nw(t + 2). This is a contradiction 
with our assumption u>(i + 1) rx co(i + 2) which means a >(i + 1) H u(i + 2) ^ 0. 
Therefore, if 9J1 is not tetrahedral, the case u{i - 1) f~M'i + 1) = w(t + 1) Dcu(t + 2) 
is not possible. 7 

^Obviously, Chore are proper self-avoiding walks for tetrahedral meshes. Note that a tetrahedral mesh cannot 
be a proper submesh of a larger 2D mesh which is a manifold as a topological space. 
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Case II: fw(i — 1) H w(i) rv o>(* + 1) 

Figure 11 illustrates the modifications necessary if wsm- T enters t through an edge and leaves 
it through a vertex. 





Figure 11: Existing walk enters the triangle adjacent to r through an edge and leaves it through 
a vertex. 

• If a >{i) n a >(i + 1) C w(i) f~l r, then a >(i) hr u>{i + 1) is a proper extension of u> (cf. 
Fig. 11(a)). 

. Suppose W (i)nw(i + 1) ^(i)flr. (Note that «(t)n«(* + l) definitely consists of a single 
vertex!) 

- If L0(i - 1) Dr ^ U)(i - 2) n u(i - 1), then w(i - 1) -> t h lo(i) is a proper extension 
of W (cf. Fig. 11(b)). 

- If w(t - 1) n r = w(i - 2) n u(i - 1), then u(i -2) a«(»- 1). In that case, 
w (t _ 2)^rh a >(*) h w(t - 1) -> w(i + 1) is a proper extension of u ; (cf. Fig. 11(c)). 


Case III: u>(i — 1) f~\ oi(i) h u>(i +1) _ _ j 

Figure 12 illustrates the modifications necessary if - T enters t through a vertex and 

leaves it through an edge. 



Figure 12: Existing walk enters the triangle adjacent to r through a vertex and leaves it, through 
an edge. 
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I- to(i) is a proper extension of w (cf. 


• If u)(i - 1) n u>(i) C w(i) n r, then u{i - 1) -» r 
Fig. 12(a)). 

• Suppose u)(i — 1) H <jj(i) ^ w(i) n r. 

- If u(i + 1) n r ^ u(i + 1) fl w(i + 2), then w(i) hr -> u>{i + 1) is a proper extension 
of u) (cf. Fig. 12(b)). 

- If u>(i + 1) n r = w(i + 1) n uj{i + 2), then w(i + 1) rv w(z + 2). In that case, 

_ i) _> w (i + 1) h u>(i) hr-)w(i + 2) is a proper extension of w (cf. Fig. 12(c)). 


Case IV: oj(i - 1) r\ u>(«) rx w(i + 1) 

Figure 13 illustrates the modifications necessary if wot t enters and leaves t through a 
vertex. 


(a) 

Figure 13: Existing walk enters the triangle adjacent to r through a vertex and leaves it through 
another vertex. 




• If w(i - 1) n w(i) C w(i) n r, then w(i — 1) — ► r H w(i) is a proper extension of w (cf. 
Fig. 13(a)). 

• If w(i — 1) fl w(i) w(i) fl r, then w(f) h r — > w(i + 1) is a proper extension of w (cf. 
Fig- 13(b)). 

Since we have been able to properly extend w w , T in all cases, our proposition is proved. □ 


B Proof of Proposition 4.1 

Proof: We prove Proposition 4.1 by induction over the level l(T) of T. Lemma 4.1 is the 

starting point for our induction (l(T) = 1). , „ ^ r 2 

Assume now that l(T)=n + 1. Let T' = T - (TnT n+l (^)). T' is a local refinement of a 

and l(T') = n. This technique is known as hierarchical coarsening [2] (cf. Fig. 14). 

By our induction assumption, there exists a proper self-avoiding walk w r>c which is com- 
patible with c. We show that there is a proper extension u> T ,c of w T /, c . 

T can be obtained from T' by local refinement. Let r be a triangle of T which must be 
subdivided (red or green) in order to obtain T. u?', c induces a constraint on the children o r 
or, in other words, the constraint on r is propagated to r’s children by w T ,, c . Tins situation is 
illustrated in Fig. 15. Figure 15(a) shows a local refinement T and the given constraint (aji). 
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Figure 14: Four steps of hierarchical coarsening. 


Figure 15(b) shows one step of hierarchical coarsening where we already have a proper self- 
avoiding walk lo T ',c ■ This self-avoiding walk induces constraints on the higher-level triangles 

shown in Fig. 15(c)). 




Figure 15: Two steps of constraint propagation. 


Since r will be refined red or green, we can apply Lemma 4.1 (and the remark thereafter) 

to properly extend lot\ c over the children of r. . , , f ™ 

If we apply this procedure to all triangles t G T which must be subdivided to get to 2 , we 

obtain u>r jC . 

This proves our proposition. □ 

Figure 16 shows an example for the correct extension of a proper self-avoiding walk to a 
higher level of refinement. Figure 16(a) shows the mechanical way dictated by the constrain s^ 
This walk can often be improved in the sense that certain jumps over vertices can be replaced 
by jumps over edges since the triangles under consideration share an edge (cf. Fig. ( ))• 




Figure 16: Extension over a higher level of refinement. 


